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Abstract 

We find a one-parameter family of long-lived physical string states 
in type II superstring theory. We compute the decay rate by an exact 
numerical evaluation of the imaginary part of the one-loop propagator. 
Remarkably, the lifetime rapidly increases with the mass. We find a 
power-law dependence of the form T = const. g~ 2 (mass) Q , where the 
value of a depends on the parameter characterizing the state. For the 
most stable state in this family, one has a = 5. The dominant decay 
channel of these massive string states is by emission of soft massless 
particles. The quantum states can be viewed semiclassically as closed 
strings which cannot break during the classical evolution. 
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1 Introduction 

A basic feature of Superstring theory is the presence of an infinite tower 
of massive states, with an exponential degeneracy at high mass level. All 
these states are expected to be unstable, since they can decay by emission 
of massive or light particles. An important question is what is the lifetime 
of massive states in superstring theory. This question was investigated in a 
number of papers in the past, using different methods P - Neverthe- 
less, there are no conclusive results and the whole picture is far from clear. 
Could there be long-lived states in string theory? The possible existence of 
quasi stable states in string theory is an exciting subject, since it could have 
both theoretical as well as phenomenological implications in Cosmology or 
accelerator physics in models with large extra dimensions. 

How to identify possible long-lived string states among the vast number 
of states with a given mass? 

In section 2 we will give a simple semiclassical argument which allows to 
find canditates for quasi stable quantum states. The basic idea of the argu- 
ment is as follows. Highly excited string states have large quantum numbers 
and can be semiclassically described in terms of solitons. Such solutions are 
classically stable, unless the string breaks. However, for a theory of closed 
string only, like type II superstrings, classical breaking is possible only if 
two points of the string touch during the evolution. Thus the problem of 
identifying stable quantum states is mapped to the problem of finding classi- 
cal string solutions whose points do not touch during the evolution. Given a 
classical solution satisfying this requirement, one can write the corresponding 
quantum state, and compute the lifetime explicitly by the quantum one-loop 
calculation. When the classical breaking is not possible, this does not pre- 
vent that the original string can decay by emitting light particles (actually 
massless, as we will see), for which a semiclassical description does not apply. 
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With this line of argument, we have identified a one-parameter family of 
long-lived string states. 

The method used here is a direct numerical evaluation of the formula 
giving the imaginary part of AAf 2 at one loop. The method is in fact similar 
to the one used by Okada and Tsuchiya in jl] for the open string state 
with maximum angular momentum. In that case, it was found a lifetime 
T = const. M _1 , implying that the open string becomes more unstable for 
large masses. We have reproduced this result as a check with our method. 
However, the states we consider here have no analog in open superstring 
theory, except for the less stable states of the family which correspond to 
classically breakable configurations. 

In sect. 3 we explain the computation of Im(AM 2 ), which leads us to a 
formula involving a finite sum which is to be evaluated numerically. In this 
paper we consider type II superstrings in 10 flat uncompactified dimensions. 

In section 4 we discuss the results of the numerical analysis. 

In the case of those states whose classical analog is a breakable string we 
find that the lifetime can increase or decrease with increasing mass, depend- 
ing on the actual quantum realization; the decay channel into two massive 
pieces is dominant and the numerical results accurately reproduce what is 
expected from the classical picture of the string breaking. 

In the other cases (of quantum states which classically correspond to 
unbreakable strings) the by far dominant decay channel is the emission of 
massless particles, mostly of very low energy. We find that the lifetime 
grows as a power of the mass, and we compute the spectrum of the emitted 
massless particles. 

In section 5 we make some additional comments, in particular we discuss 
possible gravitational effects, which were not taken in account before, and 
argue that they are not expected to play an important role. 

Some details of the calculation are described in appendices. In appendix 
A, we describe the construction of the quantum superstring states whose 
decay properties are studied in this paper. In appendix B, we construct the 
corresponding vertex operator. In appendix C we compute the correlators 
appearing in the two-point amplitudes in terms of theta-functions of the torus 
and organize the result in a convenient way. In appendix D we compute the 
fermion correlators and perform the sum over spin structures. In appendix 
E we give a simple analytic computation of a decay in a particular channel, 
obtaining a decay rate which agrees with the numerical results. In appendix 
F we show the shape of the logarithm of the decay rate divided by 2M 2 for 
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different values of M for some of the analysed cases. 

2 String solitons and quantum states 

Quantum states with high quantum numbers can be described in terms of 
a classical soliton. Examples on how a classical description can describe 
features of the quantum decay with high accuracy were given in ref. (see 
also sect. 3). In this section we construct classical soliton solutions which can 
be associated with the quantum states whose decay rates will be computed 
in the next sections. A classical soliton description of a quantum state can 
give an intuition on the possible decay modes of the quantum state. Here we 
will consider several examples and identify some quantum states which are 
expected to be long-lived. 

2.1 Quantum Hilbert space 

We review here some standard facts of the free string in ten flat spacetime 
dimensions and set the notation (covariant super-vertices will be discussed 
in the next section and in the appendices A,B). Define U = —Xq + Xg, 
V = X + X 9 , and set U = 1\fa i 'p v T. Consider strings which, in addition of 
moving in the X 9 direction, fluctuate in the planes 



We shall consider states in the NS-NS sector and omit in the whole section 2 
the world-sheet fermions from the formulas (their inclusion is straightforward, 
but not very relevant for the discussion of this section). 
The solution to the equations of motion is 



V = V L {a + ) + V R {a-) , X\a,r)=Xl(a + )+X^a-) , z = l,...,4, (1) 



where <r± = a ± r. In addition, one has to solve the two constraints: 



Z 



i 



X 1 + iX 2 
V2 



z 2 



X 3 + iX 4 
V2 
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We will work in the center of mass frame of the X 1 where Pi = 0. For the 
complex coordinates, the Fourier mode expansions is (cr G [0, 7r) ) 



g-2in(o--r)j 



n=l * 



, oo 



2 ^ x/n 1 

n=l 
, oo 

a' \ -v 1 



2in(a-~r) _ f -2in.(o--r)] 



^ = yyE^^- e " 2m(ff+T) -^ 2m(CT+T) ] (4) 

n=l * 

The operators b, c obey the usual commutation relations: 

[bn±, ^m±] = 3nm i \bn±, b m ±] = S nm , (5) 
[Cn±; ^mi] Sum j [Cn±; Cm±] S nm . (6) 

The mass formula is given by 

a'M 2 = 2(N R + N L )-a, N R = N L , (7) 
where the number operators are 

oo 

Nr = J2 n ( b l+ b n+ + b\jo n _ + C ] n+ C n+ + C ] n _C n _) . (8) 
n=l 

The expression for Nl is similar, with the change {b — > b, c — > c}. The 
normal ordering constant a is a = 2 in the NS-NS sector. In this section it 
will be ignored since we are interested in states with large Nr = Nl- 

The angular momentum components J i2 = J\2R + J12L1 J34 = J34R + J^l 
in the (Xi^X?) and (X3, X4) planes are 

00 00 
J12R = ^2(bl + b n+ - b ] n _b n J) , J 12L = ^2(bl + b n+ - &*_&„_) , (9) 

n=l n=l 
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and similar expressions for J34R, J34L, replacing {b — > c}. The physical 
Hilbert space is then constructed (in the light-cone gauge) as usual by apply- 
ing the creation operators to the vacuum Fock state. The operators b n+ 
raise one unit of the spin operator J 12 in the "up" direction, whereas the 
operators b n _, b n _ raise one unit of the spin operator J 12 in the "down" 
direction, and similarly for J 34 in the 34 plane. 



The classical string solutions have to satisfy the constraints A well- 
known example is the rotating string with maximum angular momentum, 
where the classical solution is given by 

Z 1 = Lcos(2a)e 2iT , U = 2^p v r , V = 2\fd Vu r , (10) 

L 2 = -2a'p u p v = a' 2 M 2 . 
The classical solution (jTUjl is obtained by setting 

2v a 2\j a 

and bf_ = 6^— = 5 bf- = ^1- = • With these classical values for the 
Fourier coefficients we get 

M=- t , J=f^, (12) 
a' 2a' 

which are the correct values of the mass M and spin J for the state of 
maximum angular momentum. 

To have the same values of Jr, Jl, the corresponding quantum state must 
be of the form 

1*0 = (b\ + ) N (b\ + ) N \0) . (13) 

This has 

N R = N L = N . 
Since, classically, Nr — j^j, we have to set 

L 2 

The soliton description applies in the large L limit and ^7 does not need 
to be an integer. In general the soliton is aproximately described by the 
quantum state (|13j) . with N being the closest integer. 
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2.2 Quantum string states and solitons fluctuating in 
two planes 

In what follows we only make use of the creation operators b\ + , c\ + and b\ + , 
c\ + , so for clarity in the notation we define W = b\ + , etc. (i.e. we omit the 
subindex {1+}). 

In this paper we shall be interested in the following one-parameter family 
of quantum states: 

l*M> = {k + n) \ k . n)l (6 t ) fc+ "(c t ) fc -"(6 t ) fc -"(5 t ) fc+ "|0) , <$ n |$,„> = 1 , 

(14) 

It has 

N R = N L = 2k , a'M 2 = 8k , (15) 

Ji2R = k + n, JnL = k-n, J UR = k-n, J UL = k + n, (16) 

so that J\2 = J34 = 2k (the full superstring state including world-sheet 
fermion modes is constructed in appendix A). 

The soliton solution with the same values of Jur, Jiil, J34R, Jul is as 
follows (here a' = 1) 

Z 1R = y/kT^ e~ 2i{a - T) , Z 1L = Vk^l e 2l{(J+r) , (17) 

Z m = Vk^i e - 2l{ °- r) , Z 2L = VkT^ e 2l(CT+r) , (18) 

or 



Z 1 = e 2iT (VkT^ e~ 2ia + Vk^l e 2ia ) , 

Z 2 = e 2iT (Vk^ e~ 2ia + VkT^ e 2ia ) . (19) 

It describes a rotating ellipse. It is useful to view it in a rotated frame 

Z[ = Zl + Z2 = ^2L ie 2lT cos(2<r) , (20) 
v2 

Z ' 2 = ~ Z i + Z 2 = y^L 2 e 2iT sin (2a) , (21) 
zy2 



Li = Vk + n + \/k — n , L 2 = Vk + n — y/k — n 
Two special cases are n = and n = k. 
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For n = k, one has L\ = L 2 and the classical configuration is a string 
whose projection in both the Z\^ 2 planes is a circle. This circular string 
rotates around its center like a wheel, so that the classical distribution of 
string matter is stationary. 

For n = the solution becomes a straight rotating string. Indeed, the 
resulting solution 

n = : Z 1 = Z 2 = 2\fke 2iT cos(2a) , (22) 

is classically equivalent to the string rotating in one-plane, as is clear in the 
rotated frame Z[ = 2V2ke 2iT cos(2a), Z 2 = (cf. eq. ifTTJjl ). 
However, at quantum level, the state 

l<M = T^p (6t)*( C t)*(6t)*(5t)*| > , (23) 
is physically inequivalent from the state 

l*J~> = ^yy (b'T(Py k \0) , (24) 

representing a string rotating in a single plane, with Nr = Nl = 2k and 
maximal angular momentum, which we call J ma x state. 

Indeed, the state (jUj) is an eigenstate of J 2 = ^2i< i< j< 9 (JijR + JijL) 2 
with eigenvalue J 2 = 16k 2 + 28k (maximal angular momentum for the given 
mass). The states (jl4j) however are not eigenstates of J 2 . In fact 

J 2 |$ M ) = (12fc 2 -4n 2 +28A;)|<l> M )+2(A;-n) 2 |$ fcin+1 )+2(fc+n) 2 |$ fcin _ 1 ) (25) 

In the large k limit, the states \$k,n), \®k,n±i) represent essentially the 
same semiclassical solution (having the same values of E, J\ 2 , J34, and 
approximately the same values of J\ 2 r, Jul and J34R, Jzal) so effectively 
J 2 ^ {12k 2 - An 2 + 28k) + 2{k - n) 2 + 2{k + nf = 16k 2 + 28k = 16k 2 , which 
agrees with the value of J 2 of the semiclassical solution. 1 

In the next section we will find that the decay rates of the state \$k,o) 
(corresponding to n — 0) share some qualitative features with the decay rates 
for the state |$^ maa: ) computed in [TT], though the quantitative details 
of the quantum decay are different. 

1 Note that all the states ()14JI and l|24(l are eigenstates of J\ and Jf with the same 
eigenvalue = Ak 2 + 14k; in the interactions, however, the conserved operator is J 2 . 
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2.3 Classically unbreakable closed strings 



An important feature of these solutions is that for n ^ the closed string 
cannot break classically in type II superstring theory. This splitting pro- 
cess can only happen if during the evolution there is a time To where two 
points of the string get in contact, i.e. there are two values <J\,<Ji such that 
X tl (ai, To) = X^(a 2 , r ). For the straight rotating (folded) closed string (fTUJl . 
the breaking can occur at any time. In fact, it was found in jTTJ[T^] that the 
string can decay into two strings of masses Mi,M 2 , where M 2 is a function 
of Mi. This relation between Mi,M 2 which emerges in the quantum calcu- 
lation can be understood (and, in fact, accurately described) in terms of the 
semiclassical process of splitting. 

When the classical breaking is not possible, one expects that the decay 
channel into two large masses Mi,M 2 will be exponentially suppressed (one 
can interpret this by saying that breaking is possible only by tunnelling 
effect). The reason is that for large masses M, Mi, M 2 each string in the 
process should have a soliton semiclassical description, so that the classical 
approximation is expected to apply, modulo terms which are exponentially 
small in the masses. The original string can nevertheless decay by emitting 
light particles, for which a semiclassical soliton description does not apply. 

This intuition will be confirmed by the quantum calculation in the next 
section. Moreover, in the cases we consider, we will find that only the rate 
of massless particle emission is not suppressed, with a rate that decreases as 
the mass increases. The conclusion would be that the lifetime of quantum 
states associated with general classically unbreakable states should always 
be very large, since only decay into massless particles is relevant. 

It should be noted that there are numerous classical string solutions where 
the classical breaking is not possible. In particular, one can consider small 
perturbations of the rotating ellipse solutions, 



zl = ^ellipse + eOii?( ff+ ) + )), Z 2 = Z 2 e + e{z 2R (a + ) + z 2L (a )), 



given "unbreakable" classical solution one can construct the corresponding 
quantum state. For example, one may consider the solutions 



Zi = e 2tT (VkT^ e~ 2la +Vk^ e 2iu ) , Z 2 = e 2ir { e~ 2ta +VkT^ e 2lu ) , 



etc., where Z t 



-1,2 
ellipse 
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which cannot break, indicating that the corresponding quantum states with 
the same values of Jur, J\2L and J34R, J34L, 

I^W) = M (fc t ) fc+n (c t ) fc -"(6 t ) fc -"'(5 t ) fc+n '|0) , (27) 

N = [(k + n)\(k - n)\(k + n')\(k - n r )\]- 1/2 
should be long-lived. 2 

3 Computation of Im(AM 2 ) 

We consider type II Superstring Theory in ten flat uncompactified dimen- 
sions.. The expression for the one-loop mass shift for the state | with 
square mass M| n = N = 2k — 1 in units a' = 4, is (see Appendix C): 

47V 



x E (-)" mi - m2 g(fc,n;m 1 ,m 2 )(9 2 log^ 1 (2|r)) mi (a 2 log(^ 1 ( 2 |r)) m2 

mi ,1712 

(28) 

where the numerical coefficient Q(k,n;mi,rri2) has been defined in the Ap- 
pendix C, and d is an overall constant normalization coefficient, which is 
independent of k and n. It can be seen that this expression is modular 
invariant. 

In order to evaluate the imaginary part of AM| n we expand the holomor- 
phic (and antiholomorphic) factors in powers of e l2lTT and e l2lTZ (respectively 
of e~ i2wf and e~ i2wS ). 

It is convenient to define the coefficients 7(771, p, q) by the following ex- 
pansion: 

p,q 

(29) 

and similarly for the antiholomorphic coefficient / -f(N,m 2 ;p,q). 

2 The classical solution is related to (11911 by an 0(4) rotation, but the quantum states 
\&k,n.n>) and \&k,n) arc incquivalent. 



AML = e A mj d*z e-""€ 



J2nz\N 



2Tr9 1 (z\T) 
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We have to evaluate the imaginary part of the integral 
r fp T 

TT — I 21 . l2(p-p)lTT 1 -2(p+p)lTT2 

' 3+27V-mi-m 2 
'2 

x / d 2 z e - 4N ^- 2 (v-p+Q-p- 2N ) n y+ i2 (<i~p-Q+p) 7TX (-39) 



The integrations over T\ and over x set p = p and q = q. Then we use the 
formula 

( at. ,\2N-m 1 -m 2 +3/2 

= c " l iVU ^ (31) 

y/NT{2N - mi - m 2 + 5/2) ' 



where c" is a numerical constant independent of N,p,q,mi,m2 and 

V iV N } K N N J ' 

is the typical phase-space function. This can be seen by comparing the 
integrand in eq. (j31|) with the Schwinger parametrization of a Feynman one- 
loop diagram for a field theory with vertex ~ $ 4>i4>2, with $ representing the 
field of mass M, and 0i t 2 the fields of masses M\^- Then one sees that the 
integral leads to the one-loop correction AM 2 due to the process $ — > 0i + 02 
with 

M x 2 = g and Mf = p . 

The number 2iV — m\ — m 2 — 2 is related to the orbital angular momentum 
carried in the interaction (the M\ t2 particles having in general intrinsic spins). 

The function u = 1 — 2(^ + ^) + (|p- — j^) 2 vanishes at the border 
of phase space Mi + M 2 = M. 

In conclusion, we get for the channel q = Mf,p = M\ 

ImAMj n (p,<?) = °o 9s ^ Q(k,n-,m 1 ,m 2 )j(rn 1 ,p,q)j(rn2,p,q) 

mi,m2 

j\]-2N-m 1 -m 2 +l u j2N-m 1 -m 2 +3/2 



4 2N- mi -m 2T (2N - mi -m 2 + l) ' ^ 

where cq is a numerical constant independent of k,n,p,q, which we conven- 
tionally take as c = 32(2vr) 3 . We also recall M 2 = N = 2k - 1 (a = 4 
units). 
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The decay rate for a given channel into particles of masses Mi and M 2 is 
given by 

lmAM% n (p,q) 
RkAPi i) = — ^7= — • ( 33 ) 



This includes the contributions from all final states with the same masses 
Mi,M 2 . 

The total decay rate is -R^ al = ^ p q Rk, n (p, q) and the lifetime of the 
Superstring state \$k,n) is 

T M = • (34) 

A final comment. Both Q(k, n; mi, m 2 ) and 7(771, p, q) are integers which 
become very large in absolute value by increasing k. They have been dealt 
with by means of computer programs which manipulate integers without 
approximations. In this sense we have performed an exact computation. 
The limitation comes from the machine size and we stopped our computation 
at mass M 2 = 129 for the most stable state and M 2 = 99 for the others. 
We will see that this appears to be enough for determining the asymptotic 
(power-law) formula of the lifetime. 



4 Decay rates and lifetimes 

We report here the main results for the decay M — > Mi + M 2 . 
(Conventions: c = 32(27r) 3 , a' = 4 ). 



Table 1: Decay rates: Sum over massive channels Mi j2 ^ 



M 2 


| 'J max \ 


l$fc,o) 






19 


3.450 


2.353 


0.373 


5* 10~ 5 


39 


4.713 


2.317 


0.097 


1.19 * 10~ 8 


59 


5.321 


2.183 


0.033 


2.26* 10~ n 


79 


5.718 


2.070 


0.014 


1.15* 10~ 13 


99 


6.013 


1.978 


0.007 


1.08* 10~ 15 


119 








1.6 * 10~ 17 


129 








2.2 * 10" 1S 
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Table 2: Decay rates: Sum over radiation channels Mi or M2 = 



M 2 


| (|) J max \ 


l$k,o> 






19 


7.83 


7.02 


2.91 


0.249 


39 


6.19 


4.31 


1.27 


0.050 


59 


5.46 


3.32 


0.77 


0.019 


79 


5.02 


2.78 


0.53 


0.010 


99 


4.70 


2.44 


0.39 


0.0056 


119 








0.0036 


129 








0.003 



Table 3: Lifetimes 



M 2 


| (^) <J max \ 


l$fc,o> 






19 


0.0886 


0.107 


0.304 


4.010 


39 


0.0917 


0.151 


0.729 


19.830 


59 


0.0927 


0.181 


1.246 


51.912 


79 


0.0931 


0.206 


1.840 


103.711 


89 








137.934 


99 


0.0933 


0.227 


2.510 


178.155 


109 








224.689 


119 








277, 834 


129 








337.877 
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4.1 Case n^O 

We have considered explicitly the maximal n = k case and the case n = k/2 
(the latter is expected to illustrate a generic case < n < k, with n/k not 
n/k < 1). 

In agreement with the fact that in this case the string does not break clas- 
sically, the rate in the channels with both M 1; M 2 7^ (massive channels) 
is much suppressed (for n = k, the "maximally unbreakable" string, it is 
completely negligible) compared to the channel in which Mi or M2 is equal 
to zero (radiation channels), see Tables [0 and EJ In particular, from Table 
H one observes that the sum of the decay rates for all the channels with both 
M li2 7^ goes very rapidly (exponentially) to zero. 

Therefore the dominant decay mode is by emission of a massless particle 
(Table©. 

Figure 1 shows the decay rate for Mi = as a function of Mf , in the 
case n = k = 40. We see that it is strongly peaked for M| = M 2 — / with I 
finite and small with respect to M 2 = 2k — 1. Note that this figure gives the 
spectrum of the emitted massless particle, since its energy is 



2M 

From this we see that the emitted massless particle is soft. 

As for the other decay product, the massive one, we think that it is 
likely that it is classically similar to the decaying state, although quantum 
mechanically it can well be different (see also Appendix lEj) . 

0.002r 

0.0015- 
0.001- 
0.0005- 



20 40 60 

Figure 1: Spectrum of massless particle emission for n = k = 40, M 2 = 79. 
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(35) 



4.2 Case n = 

In this case the rates for the channels M1.2 7^ (massive channels) and for 
the channels where Mi or M2 is zero (radiation channels) are comparable. 

The massive channels Mi^ 7^ contain the case of classical breaking. In 
fact we have seen in Section 2 that n = is the limiting case where the 
elliptic classical configuration degenerates into a straight line: the classical 
closed string becomes a folded rotating string which can break at any time. 

In ref. [12], by the description of the classical process of the splitting into 
two closed strings, we obtained the following masses of the decay products: 



where a G [0, 1] is a parameter specifying the splitting point o"o = an. Eq.(jHEJ) 
defines parametrically a line Mi = Mi(M2). 

We expect that the classical breaking configuration with masses ()36|) cor- 
responds to the channel of maximum rate (for M x 2 7^ 0). Indeed the quantum 
results for the decay rates are maximal on a curve Mi = Mi(M 2 ) in the plane 
Mi, M 2 . As an illustration we plot in Figure 2 the logarithm of the decay rate 
divided by 2M 2 as a function of Mf for fixed Mf = 1, 5, 10, 20 at M 2 = 99. 

Figure 3 displays the line of maximum value of the decay rate (in the 
plane Mi/M, M2/M), showing agreement with the classical curve (j3Tlj) . This 
can be compared with the similar curve obtained in using saddle-point 
method, which applies at very large M (see also |12j). 

Near the end points of the classical curve (Mi = or M 2 = 0) the process 
of light or massless particles emission is not described by the classical string 
breaking. Indeed we see in Figure 2 that the rates for M| = 1 are distributed 
in a rather flat way and it is difficult to see a definite maximum. 

The total rate for the sum of all the massive channels where both Mi 2 7^ 
is seen to decrease with increasing M (Table [H case \$k,o)), but less fastly 
than the sum of the rates of the radiation channels with Mi = or M2 = 
(TableEl case \&k,o))- This is found quantitatively in Sect. 4. 3 by a fit of the 
data with a power like formula. 

Another interesting question is to investigate, among the different mas- 
sive channels with both Mi and M 2 different from zero, which ones gives 
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Figure 2: M 2 = 99, logarithm of the decay rate divided by 2M 2 as a function of M\. 



for Mf 



1 (flattest), 5, 10, 20 (increasingly narrower) . 



the dominant contribution to the decay rate. Figure 2 suggests that light- 
particle emission could be dominant over other massive channels, and that 
the behavior for the sum of all massive channels is likely to result from a 
complex interplay of different individual behaviors, the phase space features 
playing an important role. 

For completeness we have reported in the Tables also the results for the 
Jmax quantum state. 

The three-dimensional plot of So - the logarithm of the decay rate divided 
by 2M 2 - as a function of Mi, M2 for the cases n = and J max looks qualita- 
tively the same as the one for the J max state reported in Figure 4 of [TT] . One 
interesting question is if the maxima of the curve of Figure 3 approach So = 
at large M. This is the expected result, since in the n = and J max cases 
the semiclassical breaking process should not be exponentially suppressed. 
In a more detailed view we see that So is rather smaller than zero for small 
M and shows a number of structures (see Appendix |FJ): a peak for M\ = M2 
and shallow maxima towards light values of M\^- By increasing M, those 
structures tend to disappear. In the case J m ax the region of the maxima of 
5*0 get closer to zero, suggesting that So max — > for M — > 00. In the case 
n — 0, this is less evident. 

In refs. [TT] and [12] the study of the J ma x state was based on a saddle- 
point/WKB-type approximation. This gave an accurate determination of 
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Figure 3: M 2 = 99. The continuum line corresponds to the classical Mi = M 1 (M 2 ) 
breaking curve l|3tj|) . The dots are the values of Mi , M2 for which the rate of the quantum 



the exponential part of the formula for the decay rate, giving S 0max = in 
the M — > 00 limit (precisely for the classical Mi = Mi(M 2 ) curve (|36|) ) but 
a less precise indication of the power behavior of the prefactor, and a = 1 
was reported for the lifetime T ~ M a . We will see in Section 4.3 that the 
precise power behavior of the lifetime as determined by the exact quantum 
calculation is different, and that there is indeed a long-lived rotating string 
state |3>fc,o), although different at quantum level from the J ma x state. 



4.3 Analysis of the results 

To analize in a more detailed way the dominant contribution to the total 
decay rate, it is convenient to study the different channels separately. The 
reason is that the decay rate in the channel in which a massless particle is 
emitted has a different behavior as compared to the channel in which two 

4 Since the maximum rate for Mi = M2 occurs at Mf = 14, we plot the points where 
M2 gives the maximum for < M\ < 14 and viceversa, for a total of 29 points. 



decay is maximum. 4 
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massive particles are emitted. Considering the cases separately allows for a 
more accurate convergence to the correct trend for M — > oo. [This separation 
is relevant only in the n = and J max cases, where the decay rate to two 
massive particles is significant]. 

We fit the decay rates with a power-law behaviour, except for the massive 
channels in the cases \<&k,k/2)> \®k,k) which go very fast to zero. 

To account for subleading corrections, we have fitted the decay rates 
assuming a dependence on the mass M as follows: 

/W(M) = c M~ a eM « cM~ a (l + 0/M) (37) 

or 

/ r (2) (M) = c M~ a e^ « cM~ a (l + f3/M 2 ). (38) 

We considered a favoured fit the one with minimum square deviation. The 
analysis has been carried through considering the logarithms of the data, and 
we take as the minimum square deviation v 2 = '^2 i (log[r(M i )]—log[f r (M i )]) 2 , 
where the appropriate r(Mj) can be read from Tables 




Figure 4: The dots represent the logarithm of the lifetime in the n — k ("circular") case 
vs. the logarithm of M 2 , for the different values of Table 3. The continuous line is based 
on the fit logT = ^ log M 2 + - 6.35 (with v 2 = 4 • 10" 10 ), sec also |g3} . 

a) In the J m ax case (rotating string in one plane, maximal angular mo- 
mentum) the decay is dominated by the sum of massive channels, containing 
the classical breaking. This is most favourably fitted by 

T massive = C ^ M^e^ . (39) 



19 



with v 2 = 2 • 10~ 6 (or T m „ = c g 2 M-°- 14 e&, v 2 = 2 ■ 10~ 5 ). 

We get two comparable fits for the decay rate for the total of radiation 
channels, where at least one of the decay product is massless: 

T = c g 2 s M-°' 54 e^ or 

r = c g 2 M~ 0A6 eM (40) 

with v 2 = 8 • 10" 8 in both cases. 

b) In the case n = the sum of the massive channels (which contains the 
classical breaking) dominates the decay, with a fit 

Tmassive = C ^ M" ' 58 ^ , (41) 

with v 2 = 2 ■ 10~ 6 (or T massive = c g 2 M-°- 90 e^, v 2 = 2 ■ 10~ 5 ). 
The decay rate for the total of radiation channels is fitted by: 

r = c g 2 s M" a96 e£ (42) 

with v 2 = 7 ■ 10~ 8 (or T = c g 2 M~ 1AO e^, v 2 = 9 • 10~ 7 ). 

c) For the case n = k/2, as already mentioned, the sum of the radia- 
tion channels, in which Mi or M 2 is equal to zero, is the by far dominant 
contribution, with a fit 

To = c g 2 M" 3 - 02 e£ , (43) 

with v 2 = 4 • 10~ 5 (r = g 2 ■ c ■ M" 2 - 74 e^ gives v 2 = 7 ■ 10" 5 ). 

d) For the case n = k the sum over the radiation channels utterly pre- 
dominates (the massive channels contribution is absolutely negligible). In 
this case we have more data for the lifetime given in Table 3. Since we are 
interested in the asymptotic behavior, for the fit we only use the largest 
values of M 2 = 79, 89, 99, 109, 119, 129. We obtain 

T = c g 2 M" 4 - 98 e^ (44) 

with v 2 = 6 • 10~ 10 (r = c g 2 M~ 5 - l e~v gives v 2 = 8 ■ 10~ 10 ). 

Finally, from the above analysis, we compute the lifetime behavior with 
the mass M, that is the one given by the predominant decay channels only, 
by: 

T= 1/Y (dominant) . (45) 
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This gives the lifetime behavior for our states in the limit M — > oo. Notice 
that in the cases J max and n = a naive fit to the total lifetime given in 
Table H3 would not follow this result. In fact, the trend of the massive decay 
rate, dominating over the rest, is disguised by the fact that the sum of the 
radiation channels is of comparable magnitude. 

In conclusion our results indicate, in the large M limit, for the various 
cases, a power-like behaviour for the lifetime as: 

T = 4 M a . (46) 
9 S 

In Table 14.31 we summarize the values found for a in the various cases. The 
uncertainty on a can be estimated by comparing different fits, as it has been 
reported above. For the most interesting "circular" case n = k, the value 
of a = 4.98 has a very small uncertainty, as it can also be seen from Figure 
4. Also the value of a = 3.0 for the intermediate case n = k/2 has a small 
uncertainty, although in this case the accuracy in the determination of a is 
lower than in the n = k case, because the numerical results for n = k/2 
are so far available only up to M 2 = 99. For the breakable cases, n = 
and Jmax, the values of a given in the Table 4.3 have a larger uncertainty, 
which is more difficult to estimate, because, as we said, there are competing 
decay channels of the same magnitude but different behavior (and, moreover, 
neither in those cases we could obtain numerical results beyond M 2 = 99). 



Cases 


a 


J max 


-0.25 


n = 


0.58 


n = k/2 


3.02 


n = k 


4.98 



Table 4: Lifetimes for the different cases: T ~ M a . See the discussion 
in the text for the uncertainty on a. 

In Appendix E we report a computation with the operator formalism 
giving the behavior with M of the decay rate for the channel: 

state{n = k) — > graviton + state{n = k — 1} (47) 

We obtain a lifetime for this mode T = c(l) M 5 for any I small with respect 
to k (that is to M 2 ). This is very near the result a = 4.98, and suggests 
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that the channels (Mf = 0) + (Mf = M 2 — I) dominate the decay (see 
Fig. 1). Note that, a priori, the power law behavior of the decay rate for 
this particular channel (jTTjl needs not be the same as the formula for the full 
decay rate, since there are many other states with Mf = M 2 — I. 

5 Discussion 

To summarize, the exact evaluation of the one-loop contribution to Im(AM 2 ) 
shows that there are massive states in type II superstring theory which 
are almost stable, with a lifetime growing with the mass as fast as T = 
const. g~ 2 M 5 . 

We have seen that these states can decay only by masless emission (other 
channels being exponentially suppressed), in agreement with the interpreta- 
tion that they cannot break classically. We repeat that these are superstring 
type II states. Indeed they do not exist in type I superstring theory, due to 
the fact that the type I superstring is unoriented. 

A natural question is how the decay rates computed here are corrected 
by higher loops. One particular sector of higher loop corrections are grav- 
itational effects, which are expected to be of order g 2 M. Therefore they 
should not affect the formula for the lifetime T = const. g~ 2 M 5 as long as 
M <C (in units of a'). This indicates that the lifetime can be as large 
as T = const. g~ 12 ! 

For states with masses larger than 1/g 2 , gravitational effects can no longer 
be ignored. However, note that none of these states can become black holes 
at larger masses. The reason is that the spatial extension of these states 
grows linearly with the mass, i.e. L ~ M/a! (see sect. 2), whereas in ten 
dimensions the Schwarzschild radius grows with the mass as i? sc h ~ M 1 / 7 , 
i.e. for large masses it is always the case that the size L is much larger than 
i? sc h. As M is increased, the gravitational field near a segment of the string 
becomes strong, and the question is what happens to the string configuration 
at larger couplings. In particular, the most stable state of the family, the 
n = k state representing a rotating circular string, seems a very robust string 
configuration and it is plausible that it may survive in the strong coupling 
regime. This state (or a mixture of states associated with small fluctuations 
f!26|) near this state) is a natural candidate for becoming a black ring ^3] m 
the strong coupling regime. The fact that this rotating n = k string solution 
exists in D = 5 but not in D = 4 is consistent with the fact that there 
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is no black ring in four dimensions. It would be interesting to study the 
correspondence principle for this case. 

A question of interest is how many states -among the exponential number 
of states existing at each mass level- are long-lived. This is obviously a 
difficult question to answer in general, without an explicit computation case 
by case. Nevertheless, one can try to estimate the number of classically 
unbreakable states, since for such states the only relevant decay channel is 
by emission of light particles. Like in the examples discussed in this paper, 
they are expected to have a longer lifetime. 

It should be noted that the higher is the dimension, the less is probable 
that two points of the string get in contact (e.g. viewing the average string 
state as a random walk process). This suggests that in ten dimensions, a type 
II superstring state has little chance to break. Therefore, for most states 
the decay would dominantly occur through energy leakage by radiation of 
massless (or light) modes (rather than decay into two very massive particles, 
which requires breaking of the string). 

We have found that in the classically breakable n = case the lifetime 
grows as T ~ g~ 2 M a58 , and that this is determined by the asymptotically 
dominant channel of string breaking. This would suggest that a generic 
closed string state should also become more stable for larger masses. Indeed, 
in the n = case, the closed string is folded with all points in contact at all 
times, so it can break at any time. Instead, a generic closed string solution 
can only break at a discrete set of specific times and at the specific place where 
two points touch. For this reason, the decay rate into two massive particles 
of a generic state should be lower than in the n = case. Of course, it is 
possible that for some states the decay channel by massless particle emission 
is more relevant, with the effect of reducing the lifetime as compared to the 
n = case. 

It seems unlikely that the family of states ffT^j) we have considered here 
constitutes the only states that survive the large mass limit (in particular, we 
expect that the states of the larger family (J27|) are also long-lived). Long-lived 
states are very visible among the states produced in a high energy collision, 
and it is plausible that there is a vast sector of long-lived states in type II 
superstring theory which could be produced, i.e. a vast sector of states whose 
lifetime increases with the mass. This is an exciting prospect which deserves 
further investigation. 
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A Quantum state construction 

As specified in section 2 we are interested in one particular family of super- 
string states. In eq. (J 14)) we introduced only the bosonic part of our states. 
Here we construct the full Superstring state. 

Define Z l and Z 2 as in section 2, with the expansions 0, ©• The 
bosonic contribution to the Superstring state is as in eq. (|14j). which we 
repeat for convenience: 

\®k,n) = a n (6+) fe+ "(c+) fe - n (6t)^(ct) fe +«|0) (48) 

where normalization factor: 

an = (k + n )\{k-n)\ (49) 

Now we proceed to the construction of the superstate. We look for a state 
with the same values of energy and angular momentum components. The 
superstate satisfying these requirements is given by: 

\$k,n)=Ma n \<f>«M) (50) 

with 

(51) 

where we have written only the right-moving part (the left-moving is obtained 
by substituting all the oscillators with tilded operators and changing n —* 
—n). It is obtained from (}4*%j) by applying the mode 

Cl/2 =^<% |-n 
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of the superconformal current (the analog must be done in the left-moving 
sector). Here we have used real coordinate modes in order to simplify the 
notation. 

The state (j51|) satisfies all the physical constraints, namely (we write only 
those of the right-moving sector) 

GV|$fc,n) =0, r>l/2, 

L n \®k,n) =0, U>1, 

(Lo-i)|$ M ) =0, (52) 

where L n , G r are the super- Virasoro algebra generators and we are in the 
Neveu-Schwarz sector (the zero-point energy is equal to |). 

The normalization constant in front of the state can be easily computed 
by requiring (<&*.,„ | = 1, giving: 

M=±, (53) 



i.e. an additional factor ^ with respect to the bosonic state 

Finally, one can check that the mass-squared operator M 2 acting on this 
state gives 

a'M 2 = 4(2fc - 1) . (54) 



B Vertex operators 

Here we construct the vertex operators corresponding to the state we have 
introduced in the preceding section. From now on we will consider explic- 
itly only the right-moving (holomorphic) part of the complete vertex. The 
construction of the left-moving (anti-holomorphic) part is similar. 

The vertex operator for a superstring state can be constructed in different 
pictures. Following the definitions of ^3] , we shall use the "integrated" form 
of the vertex operator. 

Given an operator W of conformal dimension hyy = |, a vertex operator 
V of conformal dimension h = 1 can be obtained by writing 

V = G_i ■ W , (55) 
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where U A ■ B " indicates the usual commutator (or anticommutator) of the 
operators A,B. 5 

As operator W we take 

W k>n = ((k+n)ip zl dZ 2 + (k-n)ifj Z2 dZ 1 )dZ 1 k+n ~ 1 dZ 2 k ~ n ~ 1 (rightmovers)e ip ' x , 

(56) 

which has conformal dimension hw = \ provided p 2 = 2k — 1. 

Acting on this operator with the super conformal charge produces a certain 
number of terms. We are interested only in the terms which will give a non- 
vanishing contribution to the one- loop amplitude. As shown in appendix ID1 
only a few terms with a specific combination of fermionic operators give a 
non-zero contribution in the two-point function on the torus. 

Considering therefore only these terms, our final vertex operator is 

V Kn = NormV L V R e ip - x , 

with 

V R = (^ Z W 2 + ^ Z2 d^ Zl )(k 2 -n 2 ){dZ 1 ) k+n - 1 {dZ 2 f- n - 1 (57) 
+ (k + n){k + n- l)i> Zl di> Zl (dZ 1 ) k+n - 2 (dZ 2 ) k - n 
+ (k - n)(k — n— V)i) Z2 d^ 2 {dZ x ) k ^ n {dZ 2 ) k ~ n - 2 . 

The Left part of the vertex V^ n is the same with d — > d and ip — > ip. 
We recall that the normalization constant is: 

4fc-4 




Norm = — — — — \ — (58) 

2k(k + n)\(k - n)\ 1 



In what follows we choose units where a 1 = 4. 



C Computation of the amplitude 

The correlator (Vk, n Vk,n) on the torus is (after summing over spin structures, 
see Appendix |D|) : 

(59) 

5 An equivalent way to construct a vertex operator is to define V = G r ■ W by requiring 
V to be independent of r. We have also pursued this way of computing V finding that 
indeed it leads to the same final expression. 
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where 



Oftfc) = ^k 2 -n 2 ) 2 O R ' L (z,0) + (k + n) 2 (k + n-l) 2 Of L (z,0) 

+ (k-n) 2 {k-n-l) 2 O R ' L (z,0) , (60) 

and 

Of (z,0) = (dZ 1 ) k+n -\z)(dZ 2 ) k - n -\z)(dZ 1 ) h+n -\0)(dZ 2 ) k - n - 1 (0) , 
OjM) = (dZ 1 ) k+n -\z)(dZ 2 ) k - n (z)(dZ 1 ) k+n - 2 (0)(dZ 2 ) k - n (0) , 
O?(z,0) = (dZ 1 ) k+n (z)(dZ 2 ) k - n - 2 (z)(dZ 1 ) k+n (0)(dZ 2 ) k - n - 2 (0) , 
and the "Left" Operators L are the same as the corresponding R with 

We see that we have 9 possible contractions (0 R L ). The basic correla- 
tors are 



2 


(* 


log 














2 


^<9 2 


log 







72 



7T 



(9^(z)9Z,-(0)) = (92i(z)9Z,-(0)) = <%2 - , 

/win 4(2fc-l) 



-ipX(z) 



,tpX(0)\ 



-4(2fc-l)^ 



6»i(0h 



(61) 



(here y = lm(z)). The generic term (0 R L ) gives (with the appropriate 

m, n, in, n): 

{ (dz,y"(z)(dz,)"(0)(l>z i y n (z)(az,f(id) > 



m!n!m!n! 



j\j\(m - j)\(m 



{dz^zmy^H^zmy 



x (dZiWdZ^O))™-* (82^)3^(0))™-* (62) 

Note that in every term the relations m — j = n — j , m — j = n — j hold, 
and thus j = j — (m — n) = j — (n — fh). We finally obtain: 

4(2fc-l) 

(v fc) „(z)y M (o)) 



(A; + n)! 2 (A;-r2)! 



4A; 2 



^i(o|t; 

J-0-2 



-4(2fc-l)7ri; z 

e T 2 



<9 2 log^^|r) + ^ 

7T 



2(j+0 



27 



The sum runs over 2n — 2 < j < &; + n, < I < k — n, 

TA=iSi{j,l]k : n) 



c(j,l; k,n) 
and 



j\l\(k + n- j)\ 2 (k -n- l)\ 2 (j -2n + 2)!(Z + 2n + 2)! 

3 

^2 ^C?' hk,n) = (k + n- j) 2 (k -n- l) 2 A(j, /; k, n) , 
i=i 

6 

^ SiO', Z; A;, n) = (k + n - j) 2 (k + n-j- l) 2 B(j, I; k, n) , 

9 6 

^Si(j,l; k,n) = ^Si(l,j; k, -n) , 



(63) 



i=4 



i=7 



i=4 



with 



A(j,l;k,n) = A(j - 2n + 2)(Z + 2n + 2)(Z + 2n + 1)(Z + 2n) 
+ 4(j - 2n + 2)(j - 2n + l)(j - 2n)(Z + 2n + 2) 
+ 16(j - 2n + 2)(j - 2n + 1)(Z + 2n + 2)(Z + 2n + 1) ,(64) 

B(j,l;k,n) = 4(j - 2ra + 2)(Z + 2n + 2)(Z + 2n + 1)(Z + 2n) 
+ (Z + 2n + 2)(Z + 2n + l)(Z + 2r2)(/ + 2n- 1) 
+ (j -2n + 2)(j - 2n + 1)(Z + 2n + 2)(Z + 2n + 1) . (65) 

Note that for the particular case n = k the above expressions simplify, 
since the only nonvanishing contribution occurs for j = 2k — 2, / = and 
c(2k - 2, 0; k, k) = l/(2fc - 2)! 2 . 

Expanding the holomorphic and antiholomorphic binomial factors we ob- 
tain 

4(2fc-l) 



(V k , n (z)V k , n (0)) = e 



-4(2fc-l)^ 



0i(z|t) 



6»i(0|r) 



7T 



r 2 



2(2fc-l)-2 



y (_)-m 1 -m 2 g (A . )n . mi)m2)(a 2 log ^ i ^| r)) m 1 ^2 log ^ i(z | r)) m 2 ( gg) 



mi ,r?i2 

where 



(fc + W )! 2 (fc-n)! 2 ^ 

Q{k,n;m 1 ,m 2 ) = ^2 



ij,l;k,n)(j + l)\ 2 



j.i 



mi!m2!(j + Z — mi)\(j + Z — m 2 )! 



(67) 
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In the particular case n = k we get 
Q(k, k; mi, m 2 ) 



(2k-l)\ 2 



m 1 \m 2 \(2k - 2 - mi)!(2Jfe - 2 - m 2 )! 
The computation for this n = k case can be done directly from the vertex 
1 (2k -V 2 



(68) 



k,k 



2 2k-2 (2fc-l)!2 

noticing that in this case 



2k-2 ipX 



(69) 



(V k , k (z)V k , k (0)) = (2k-l) 2 (e-^e^) 



7i 



d 2 log0i(z|r) + — 

T2 



2(2fc-2) 



• (70) 



Another state of interest is the state of maximal angular momentum with 
M 2 = 2k — 1 which rotates in one plane only. In this case the superstate is 
(cf. eq. (J24j) . for the bosonic version) 



(2k-l)\ 

The corresponding vertex is given by 

Jmax 1 — 1^ 



^i + (6 t ) 2ife - 1 ^i(6 t ) 2fc - 1 |0> . 



and 



2 2fc-2 (2A; - 1)! 



(V k Jmax (z)V k Jmax (0)) = ( e ^^V pA(U) ) 



2fc-2 ipX 



(71) 



(72) 



ipX(z ipX(0)\ 



X 



(2A; — 



i! 2( 2fc _ 2 _ i)! 2 



^ log0i(«|r) 



7T 



t 2 



2/ 



7T 



'r 2 ' 



2(2fc-l)-2j-2 



One obtains 



g J " iaa: (fc ) m 1 ,m 2 ) 



(2Jfe - 1)! 2 (2(2£; -I) - my- m 2 )! 
m 1 \m 2 \(2k - 2 - mi)! 2 (2£; - 2 - m 2 )! 2 



(73) 



(74) 



In all cases AM 2 is finally obtained by the integral: 



AM 2 



S / d 2 , (V(^)V(O)) 



(75) 
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by inserting the appropriate vertex. 

We will use convenient formulas for the expression involved in the com- 
putation of the amplitude IT5|) 

2 ^'^W) = (76) 



\ni3 1 t\n-l„n(n-l)/2 (n-l)(n-2)/2 , 2 n-l 



v.n=l / \n=l 



oo oo 

— ^(log^r)) = ^n^+^M n, (^ + ^) (77) 

n=l n,i=l 

with xi = e 27ri(T " z) , x 2 = e 2 ™ 2 . 



D Sum over spin structures 

In the NS-R formulation of type II superstring theory, the product of func- 
tional determinants of bosonic, fermionic coordinates and ghosts for a definite 
spin structure s is (Right sector) 



The other ingredient which depends on the spin strucure s is the fermionic 
correlator: 

Summing over spin structures with the GSO phases rj s , we get 
Since #i(0) = we obtain a vanishing result for z — > w. 
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Now consider correlators involving dip. Making first the derivative in z 
and w we get a non- vanishing result for z — ■> w: 

linw ^^Z s (^(z)^(0)) s (* A H^ p (0)) s 

s 

9 " v9xp {B 1 {z)f{9'Mf~ 9 '" v9Xp - [ ] 

This is the basic correlator used in appendix in order to obtain eq. (|59j). 



E Analytic determination of a decay rate 

We study here the decay of the state n = k in a particular channel where 
one of the decay products is massless, by using the operator formalism. 

Let us for short use the following notation for the state n = k (in this 
notation N = 2k — 1): 

\$n) = ^ 4(b j ) N 4(c j ) N \0) ■ (80) 

We will compute the decay-rate of the process in which the state n = k 
decays into a massless particle ( "graviton" ) and another massive state of the 
same kind: 

|$jv) -> \graviton) + \$ N -i) (81) 

We are interested in the behavior of the decay rate for N ^> 1 and I finite. 
Thus we consider the matrix element 

A = g s A R A L (82) 

where 

A R =(0\-^===ij 1 e-bexp(-V2p-b) ^=|0) (83) 

Here e is a polarization tensor and in complex coordinate notation e-b = 
e\b and similarly p ■ b = p\b. 

We have written only the part of the graviton vertex which is most rel- 
evant for the N — >• oo limit. Other parts will be suppressed by additional 
powers of p\ < p, the graviton energy being 

I 

P = 7= ■ 
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The computation of the matrix element gives: 



Ar - €i (/-i)! \jw^iy-' { ) 

and similarly for Al with ei,pi — > 62,^2- 

In the limit N — > 00 at fixed /, — > c(/, #) v^/V where c(l, 9) depends 
on the angle of the graviton momentum with the (complex) directions 1,2 
respectively. 

We thus get the decay rate for the present channel: 

D(N, I) = -±=g 2 s J cM 9 \A R ■ , (85) 

where the last phase space factor comes from 

f , ^ 6(VN-p-y/^ + N-l) = ^=, 

J Pv /p 2 + N-l s/N 

and we have ignored some numerical factors. The integration over the solid 
angle gives a suppression factor which depends on I but not on N. 
Finally we get the behavior: 

D(N, I) = g 2 s co(Z) iV" 5/2 , N » 1 . (86) 

We think that other channels are also relevant in the process M — > M\ + 
massless. In those channels the massive final state can be different from the 
class of u n = fc" states considered here. 
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F Three-dimensional plot of the logarithm of 
the decay rate for n = and J max cases 

Here we show graphics representing the shape of the logarithm of the decay 
rate (divided by 2M 2 ), called Sq in the text, as a function of Mf, M|. 

The figures have been obtained focusing on an interval in the phase space 
in which < M\ 2 < 20, and limiting the value of Sq to a range [—0.4, 0]. 

We have chosen a perspective so as to have a clear view of the distance 
between the maxima and the plane Sq = (cf. also figs. 2 and 3). 



Case n=0, M =19 Case Jmax, M =19 




Figure 5: Three-dimensional plot of Sq. 
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